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Abstract 

Spin-torque oscillators are strong candidates as nano-scale microwave generators and detectors. 
However, because of large amplitude-phase coupling (non-linearity), phase noise is enhanced over 
other linear auto-oscillators. One way to reduce nonlinearity is to use ferromagnetic layers as a 
resonator and excite them at localized spots, making a resonator-excitor pair. We investigated 
the excitation of oscillations in dipole-coupled ferromagnetic layers, driven by localized current at 
ferromagnetic nano-contacts. Oscillations possessed properties of optical-mode spin-waves and at 
low held (^200 Oe) had high frequency (15 GHz), a moderate precession amplitude (2-3°), and a 
narrow spectral linewidth (< 3 MHz) due to localized excitation at nano-contacts. Micromagnetic 
simulation showed emission of resonator’s characteristic optical-mode spin-waves from disturbances 
generated by domain-wall oscillations at nano-contacts. 

PACS numbers: 75.30.Ds,75.40.Gb,75.78.Cd,85.75.-d 
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I. INTRODUCTION 


Transfer of angular momentum from a dc spin-polarized current to a nano-scale fer- 
romagnet (FM) exerts an anti-damping torque, Spin-Transfer Torque (STT), that can 
compensate intrinsic damping torque and induce stable precession of magnetization. 1 When 
combined with a magnetoresistance effect, like giant magnetoresistance (GMR) or tunnel¬ 
ing magnetoresistance (TMR), high-frequency voltage oscillations are emitted, producing 
Spin-Torque Oscillations (STO). 2 The same structures can also rectify injected ac voltage 
at resonance) 3 ^ and off-resonance 5 . Such microwave nano-oscillators/detectors are sought 
after for applications like inter-/intra-chip communication, imaging 6 and non-destructive 
testing, and lab-on-chip sensors. STOs are suitable candidates for such a role. However, 
low output power, the trade-off between power and frequency, and weak coherence of os¬ 
cillations hinder their applications compared with their semiconductor counterparts. The 
large nonlinear phase-amplitude coupling 7 ' produces a dilemma between large precession 
amplitude and small linewidth. Also, frequency locking of STOs to external reference signal 
becomes hindered. 8 Such a dilemma can be overcome by the separation of STOs into an 
excitation source and a resonating element, so that precession frequency and amplitude will 
be set by the resonator design, not by the driving STT. Then the linewidth will be reduced 
considerably. 

It was shown that non-uniform current density in TMR-STO resulted in an increase the 
amplitude of generated precession, 9 12 and reduction in linewidth. 9 13 15 However, the origin 
is still not clear, and the fabrication process is not well understood or easily reproducible. 
On the other hand, the Ion-Assisted Oxidation (IAO) of ultra-thin aluminum reproducibly 
was used to fabricate a 1-nm-thin alumina Nano-Oxide Layer (NOL) with direct 2-nm 
Nano-Contacts (NCs) between FM layers.® ® In NOL-based STOs, moderate power and 
narrow linewidths were reported)®® 4 with oscillation behavior similar to low-TMR-STO. 25 
There is evidence for the presence of FM metallic NCs by magnetoresistance and trans¬ 
port properties,® transmission-electron micrographs) 18 1 20 1 - and conductive atomic force 
microscopy.® 19 26 However, the measured Nano-Contacts Magnetoresistance (NCMR) ratios 
are far below expectation compared with scattering from confined Domain-Walls (DW), 27 ' 29 
mostly due to the presence of oxygen and non-magnetic impurities) 19 26 3 °I 
In this paper, we propose that the localized precession of DWs at NCs work as excitors 
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of spin-waves in FM layers. This makes frequency completely determined by resonator’s 
designed eigen-frequency, regardless of the mechanism of excitors. The loss of frequency 
tunability reduces non-linearity and linewidth considerably After the experimental and 
simulation description, we present the excited modes in the chosen resonator, then we dis¬ 
cuss the reduced non-linearity of NCMR-STO. For the resonator, we used a nano-pillar with 
two free FM layers, where magnetostatic dipolar field provides inter-layer coupling with two 
coupled-oscillations characteristic modest The dynamics of coupled free FM layers are of 
practical interest for high-frequency emission at low applied field (< 500 Oe), linewidth nar¬ 
rowing, and doubling of magnetization precession frequency in resistance oscillations. 32 35 
Most importantly, the loss of current-tunability of frequency can be compensated by the 
controllability of inter-layer relative angle. 


II. EXPERIMENT AND SIMULATION 

The film stack with designed thickness in nm was: thermally-oxidized silicon sub¬ 
strate/electrode layer(Ta 5/Cu 200/Ta 40/chemical-mechanical polishing)/milling 5/Ta 
3/Ru 2/Fe5oCo5o 5/A1 1.3/IAO 20 seconds exposure time/Al 0.3/Fe5oCo5o 5/capping (Cu 
10/Ru 10). The film was deposited by magnetron and ion-beam sputtering in the chambers 
described before. 19 Subsequently, films were vacuum-annealed at 270°C and 400°C for 1.5 
hours each with 10-kOe magnetic field. The choice of capping material and additional metal 
A1 insertion were optimized with annealing process for lower Resistance-Area (RA) product 
and enhanced MR ratio, based on previous work. 19 30 Current-Perpendicular-to-Plane (CPP) 
pillars of elliptical cross-section were micro-fabricated by Ar + ion-milling and electron-beam 
lithography. RA was found from the slope of 4-probe dc resistance vs. area inverse (R-l/A) 
line, and compared with Current-In-Plane-Tunneling (CIPT) measurement of unpatterned 
films. 

The STO microwave emissions were measured by a 26-GHz spectrum analyzer under bias¬ 
ing from a bias-T, and scattering at the measurement probe was measured by a network 
analyzer. Geometry, angles and coordinates definitions are summarized in Fig. B>). The 
positive current was defined to be electrons flowing up. We are presenting the detailed 
measurements of a 320 nm x 160 nm pillar at £ = 60°, although results presented later were 
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qualitatively similar among samples. 

Micromagnetic simulation was done by Nmag, a finite-element-method simulator. 36 l Calcula¬ 
tion geometry consisted of two ellipses same as experimental design separated by 1 nm, and 
NCs were included as 1-nm-radius cylindrical contacts between the two layers. Inter-layer 
dipolar field was included through demagnetization held calculation. Material parameters 
are: stiffness constant of 2.3 x 10 -10 erg/cm 3 , saturation magnetization of 1930 emu/cm 3 
with Gilbert damping constant of 0.02, and an unphysical spin polarization of 100%. Mesh 
size away from NCs was set to 5 nm, and it changed to 0.7 nm inside NCs. For hysteresis 
loops, ^-dependence was calculated for 20 randomly-placed NCs. For qualitative understand¬ 
ing of STO dynamics, we compared the cases of zero and four NCs, under the experimental 
conditions of 250-Oe held applied at £ = 60°. The current prohle was approximated to be 
conhned in NCs with conhnement extending 1 nm away from middle of NC into FM layers, 
as most of the voltage drop will be on this region, 38 although more accurate representation 
is needed. 39 The total current was +17.5 mA and the current distribution was calculated by 
assuming that a single NC and tunnel barrier resistances are 600 fl and 500 0, respectively. 18 


III. RESULTS AND DISCUSSION 

RA found from the slope of 4-probe R-l/A line, and CIPT measurements were 0.2 and 
0.3 Cl ■ /im 2 , respectively. The bias dependence of 4-probe differential resistance at parallel 
magnetization state (inset in Fig. [IJb)) was relatively hat. The resistance temperature 
dependence of similarly conditioned hlms also showed metallic-transport character. This 
indicates that conduction is dominated by transport through NCs and not by tunneling 
through oxide barrier. Previously, high-temperature annealing (> 380°C) was hindered 
by manganese diffusion from pinning antiferromagnet towards NOL. 19 Better NOL barrier 
quality and purer NCs were obtained in this report by using manganese-free structure for 
higher annealing temperature, in addition to insertion and capping layers optimizationsJ 30 
Figure [IJb) shows the two-probe resistance vs. magnetic held (R-H) applied at £ = 0° 
and 60° measured at the same position as STO measurements. From the switching helds 
of easy-axis R-H, interlayer dipolar coupling held (H ic ) is estimated at 400 Oe. It is in 
agreement with the estimation of cross-demagnetization,® H ic = 47rpi 2 M s = 433 Oe, where 
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M s =1750 emu/cm 3 is the measured saturation magnetization, and pu = 0.0197 is the cross¬ 
demagnetization factor. This dipole inter-layer coupling can be considered equivalent to 
the usual bilinear coupling through a metal spacer (J = —2 dM s H ic = —0.7 erg/cm 2 
defined negative for anti-parallel coupling). The contribution from coupling through NCs 
and spacer roughness to magnetostatic energy can be neglected. Ferromagnetic coupling 
was found to be small from free-layer magnetization-loop shift of unpatterned spin-valve 
films (IrMn/FeCo(Pinned)/NOL-NCs/FeCo(Free)), with Jncs = 0.01 — 0.02 erg/cm 2 . 19 
Micromagnetic simulation reproduced static R-H, H ic estimation, and the reduction of AP- 
to-P plateau width by NCs (Fig. [l](c)). We chose for oscillation measurements the pillar 
that had the closest R-H curve to micromagnetic simulation, which had 11.1% MR ratio 
and 0.17 flp, m 2 RA product. Due to large pillar size, uniform rotational switching was not 
reproducible for field applied along easy axis (£ = 0°). At tilted angles, the magnetization 
rotated as a single domain. Thus we are presenting tilted angles results of oscillations 

(Fig# 

Largest power microwave oscillations were observed for £ = 60° at ~15 GHz when applying 
high currents. Sample power spectrum with a Lorentzian peak fitting is shown in Fig. UN- 
There is a drop in resistance at I<j c = 14.7 mA accompanied with a jump in oscillation 
frequency, f osc , a narrowing in full-width-at-half-maximum, A/, and increase in integrated 
power, Pi n t, indicating a change into auto-oscillation mode,®] with a mechanism similar 
to STOs based on pin-hole tunnel junctions (Fig. [2](c))P Linear fits to normalized inverse 
power, 1/p, at sub-threshold gave a threshold current, I t h of 14.74 mA. The presense of two 
frequency branches at sub-threshold and high-current regions can be ascribed to edge and 
center modes in elliptical geometries! 43 The highest oscillation power is 0.4 nW (1.6 nW if 
corrected for impedance mismatch) giving a precession amplitude (9 P ) of 2-3°, whereas the 
lowest A/ is 3 MHz corresponding to a quality factor of 5000. 

Regarding the excited oscillations, the possible coupled oscillations or spin-waves in two 
layers of free spins are the optical (anti-phase) mode (OM) and the acoustic (in-phase) 
mode (AM). 31 Modes frequencies can be found from the solution to coupled Bloch equations 
of the two layers with effective field determined from the free energy.^ 1 44 45 We considered 
only the main contributions of Zeeman energy, film demagnetization, and interlayer dipolar 
coupling. The optical and acoustic eigen-frequencies of an in-plane magnetization precession 
can be simplified to: 


5 




fa,' 


7/27T 

/op 

7/27T 


= (H cos ip + 47 tM s — 2H ic (cos A 6 + 1)) (H cos ip) + 8 Hf c cos A 6, (la) 

= ( H cos ip + 47 tM s — 2 H ic (cos A 9 + 1)) ( H cos ip — 4 H ic cos A 9 ), (lb) 


where 7/277 = 2.8 MHz/Oe is the gyromagnetic ratio, and other symbols are defined in 
Fig. [IJa). We confirmed the presence of the weaker-amplitude AM (Fig. [2jbj ). Using H tc 
= 400 Oe and A 9 = 130-150° from measured R-H and micromagnetic simulation results 
in / op of 14.0-16.1 GHz and / ac of 3.9-3.5 GHz, which agrees with the observed spectrum. 
The frequency of OM depends mostly on the coupling strength and relative angle between 


the layers (the last term on right in Eq. lb). The weak dependence of / osc against I <: ] c and 
H (~ —1.3 MHz/Oe not shown) supports that / osc is determined mainly by excitation of 
an OM spin-wave. The maximum f osc (H = 0, Ad = 180°) from other devices was 17.8 GHz 
which corresponds to Hi c = 460 Oe, in agreement with the estimation from the correspond¬ 
ing R-H curve. To increase the oscillation frequency, we reduced the size of elliptical pillars 
to 160 nm x 80 nm. At I c | c = 2.8 mA, H = 185 Oe, and £ = 50°, the resulting oscillations 
had / osc , df/dldc, and A/ of 23.3 GHz, < 4 MHz/mA, and 1.3-MHz linewidth, respectively. 
The corresponding quality factor is more than 17,000. 

Although the presented frequency of OM is higher than AM, the measured and simulated 
peak intensities of OM are much larger (Figs. [2](b) and |3|a,c)). This is due to two reasons. 
In OM, anti-phase dynamics maximize dynamic MR change. 3334 In comparison, OM am¬ 
plitude as measured by Brillouin light scattering is reduced due to canceling contributions 
to light scattering cross-section. 41 Secondly, the energy required to excite OM is smaller in 
anti-coupled (J < 0) harmonic oscillators. The average energy difference between AM and 
OM with equal precession amplitudes is: (E ac ) — (E op ) = — J (fimi • fim 2 ), where 5mi • fim 2 
is the characteristic-mode’s dimensionless power. For J < 0, excitation of AM requires 
higher energy than same-amplitude OM. 

It should be noted that Eqs. [I] were derived for infinite-wavelength limit (i.e. wave-vector- 
thickness product qd ~ 0). Quantitative corrections due to dynamic dipolar coupling 
between propagating spin-waves cannot be ignored because qd = 1.74 from simulation pre¬ 
sented later. 41 44 45 But due to experimental uncertainty in determining i7j c cos A 9, exact 
quantitative comparison becomes difficult, and the main conclusions are not changed. 

The linewidth broadening of STOs compared to linear auto-oscillators is understood to 
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be due to amplitude-phase coupling,^ which is expressed by the nonlinearity parameter 
(v)W^ The nonlinearity of presented results is v = (I&JY g )(df / dl& c ) ~ —0.16. The natu¬ 
ral FMR linewidth ( T g = 934 MHz) is obtained from linear extrapolation of A/ to zero 
current at sub-threshold J- 47 and the agility of oscillation frequency in current {d,f /d,I,\ c ) was 
-9.6 MHz/mA. This non-linearity is one order of magnitude smaller than other reported 
values? 5 47 ® Because the nonlinearity is very small, the sudden change of oscillation into a 
single mode at threshold hinders the applicability of determining v from A/'-l/p plot.® The 
Af-l/p slopes were 3.8 and 19.4 MHz/(mA 2 • //W) for above-threshold and below-threshold 
regimes, respectively. 

NCMR-STO usually showed relatively small agility (16-18 MHz/mA) 21 23 2 ® compared with 
other TMR-STOs, leading to smaller nonlinearity and narrow linewidth. Possible reasons for 
lowered agility and nonlinearity in this report can be the optimized fabrication process with 
purer NCs®® the coupled oscillations of two layers,® 35 ® and the tilted magnetization 
angle away from easy axis? 0 However, the loss of agility and small linewidths were obtained 
for various angles and frequencies, which indicates that the improved purity of NCs is the 
main factor. 

We used micromagnetic simulation to show how nano-magnets with NCs worked as a 
resonator-excitor pair. Results shown in Fig. [3] agree reasonably with experimental data. 
In the case of no NCs (upper part in Figs. [3](a— c)) oscillation frequency is similar to the 
calculated and measured ones, with the optical mode being the dominant component. But 
precession amplitude ( 6 P = 0.05°) is very small compared with the experimental value (2- 
3°). With the insertion of four NCs, optical spin-waves were emitted from NCs (Fig. [3](d) ), 
and increased 9 P to 2° (lower part in Fig. |3](a-c)). The origin of perturbation near NCs 
is of similar origin to previous reports?- 1 52 The domain-wall is pushed outside NC-region 
into the ferromagnetic layer and starts to oscillate at high frequency between Neel and 
Bloch walls (250 GHz for the chosen geometry and current density) (Fig. |3](e)). These very 
high frequency oscillations were localized up to 5-10 nm away from NC (Fig. [3](f )). The 
localized precession acts as a point source that generates spin-waves propagating radially at 
the characteristic mode of the system, which is an optical spin-wave. 

The implication on the nature of current-induced dynamics is that magnetization preces¬ 
sion is not induced by STT directly. In Fig. [2](c), at first going above I t h, localized STT 
increases amplitude, compensates damping around NCs and increases local precession am- 
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plitude. When local precession amplitude saturates, at I = 15.6 mA, increase of current 
and STT will not change 9 P , leading to loss of agility and narrowing of linewidth. At this 
stable regime, magnetic layers act as a resonator that is excited by the energy coming 
from the point sources at NCs. This makes the presented oscillator similar to a classical 
auto-oscillator, and results in a considerable reduction in linewidth. 


IV. CONCLUSION 

In conclusion, we presented the measured spin-torque-driven oscillations of a spin-torque 
oscillator with nano-contacts between two free ferromagnetic layers coupled antiferrmagneti- 
cally with dipolar-field. Resulting oscillation character agrees with propagating optical-mode 
spin-waves. In micromagnetic simulation, inclusion of NCs with localized current density 
showed that NCs work as point sources, and optical-mode spin-waves were excited in the 
ferromagnetic layers. The ferromagnetic layers act as a resonator that is decoupled from 
mechanism of excitation point sources. Since oscillation is not driven by spin-transfer torque, 
linewidth decreased. 

Although we studied dual ferromagnetic layers in this report, same effect should be the origin 
for the small non-linearity and linewidth narrowing common in NCMR and pin-hole-TMR 
STOs. 9 13 21 23 25 So, by utilizing NCs and using a magnetic resonator that has a large os¬ 
cillation amplitude, e.g. tilted-anisotropy ferromagnetsp 3 we expect to increase both power 
and quality factor for applications. 
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FIG. 1. (a) A schematic of the elliptical pillar geometry and definitions of coordinates and angles, 
(b) Two-probe resistance vs. field applied at £ = 0° and 60°. Inset shows the bias dependence of 
differential resistance in parallel magnetizations state, (c) Normalized MR found from micromag¬ 
netics simulation, for the cases of no NCs (dashed lines), and with 20 NCs (solid lines). 



13 





































FIG. 2. (a) Representative oscillation power spectrum at £ = 60° with Lorentzian peak fitting 
(solid line), (b) The acoustic and optical modes of coupled oscillations were observed near 3 GHz 
and 15 GHz, respectively, (c) Current dependencies of oscillation characteristics. Same symbols in 
/osc and A/ panels correspond to each other. 
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FIG. 3. (a-c) Comparison of micromagnetic dynamics with (top part) or without (bottom part) 
NCs at H — 250 Oe, £ = 60°, Id c = 17.5 mA . (a) Normalized magnetization’s y-component 
of the total system (m y ). (b) Magnetization angle of top and bottom layers at (x,y) = (0,0). 

(c) Spectrum of m y transformed from 250-ns (15-ns) duration of dynamics without (with) NCs. 

(d) Time snapshots with the color representing [m y (r, t) — m y (r, 0)]. Optical-mode spin-waves 
are emitted from localized excitation at NCs. (e) The origin of localized excitation is a confined 
domain-wall that is oscillating at 250-GHz. (f) Power profile of localized precession around the 
top-right NC. 


297 298 299 300 260. .1 .2 .3 




Frequency [GHz] 






Time: 0.017 ns 


5 nm 


Time: 15.622 ns 



15 



































